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' Abstract 

We exhibit a surprising relationship between elliptic gradient systems of PDEs, multi- marginal Monge- 
Kantorovich optimal transport problem, and multivariable Hardy-Littlewood inequalities. We show that 
the notion of an orientable elliptic system, conjectured in [6] to imply that (in low dimensions) solutions 
with certain monotonicity properties are essentially 1-dimensional, is equivalent to the definition of 
a compatible cost function, known to imply uniqueness and structural results for optimal measures to 
certain Monge-Kantorovich problems [11] . Orientable nonlinearities and compatible cost functions turned 
, out to be also related to submodular functions, which appear in rearrangement inequalities of Hardy- 

(-H Littlewood type. We use this equivalence to establish a decoupling result for certain solutions to elliptic 
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PDEs and show that under the orientability condition, the decoupling has additional properties, due to 
the connection to optimal transport. 

1 Introduction 

> . 

The main purpose of this note is to pinpoint a surprising connection between elliptic systems of PDEs, 
multi-marginal optimal transportation, and multivariable extended Hardy-Littlewood inequalities. A recent 
OA ! paper by Fazly and Ghoussoub [B] introduced the concept of an orientable elliptic system (Definition 12.11 

below), which seems to be the appropriate framework for investigating De Giorgi type conjectures ((?], [2]) for 
systems of more than two equations. On the other hand, the thesis of the second author [TT], following work 
of Carlier [5], introduced the concept of a compatible cost function (Definition ^. 2l below). a natural, covariant 
condition ensuring uniqueness and structural results on solutions to a multi-marginal optimal transportation 
problem with one dimensional marginals. These notions turned out to be also related to submodular (or 
2-monotone) functions, which appear in rearrangement inequalities of Hardy- Littlewood type as studied by 
several authors dating back to Lorentz [5] . 

We will show here that these three conditions are actually equivalent. As a consequence we shall see how 
the concept of an H- monotone solution to the system ([1]) below, which was introduced in [6], is intimately 
related to the geometric structure of optimal measures in the optimal transport problem @ uncovered in [llj . 
We also show that monotone solutions to the elliptic system can be decoupled. If the solution is ii-monotone, 
we use the connection with optimal transportation to show that the sum of the decoupled non-linearities is 
everywhere less than the original non-linearity. 

Let H : R m -> K be a C 2 function. We will consider the system of elliptic PDEs on R^: 

Au = VH(u), (1) 

where u — (ui,U2, ...,u m ) represents an m-tuple of functions. In this context, H is often referred to as the 
non-linearity of the system. 
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Given probability measures fj,i,fj,2, ••• > Mm on the optimal transport (or Monge-Kantorovich) problem 
consists of minimizing 

H(p 1 ,p 2 , ...,p m )dj(pi,p 2 , ...,Pm) (2) 



among probability measures 7 on W n whose 1-dimensional marginals are fa. In this setting, H is called the 
cost function. 

If H is bounded below on R m , then there exists a solution 7 to the Kantorovich problem ([2]), as well as 
an m-tuple of functions (Vi, V2, • V m ) ~ called Kantorovich potentials - such that for all i = 1, m, 

V i ( Pi ) = inf (H(pi,p i ,...,p m )-Y l Vl(Pi)), (3) 
and which maximizes the following dual problem 

m . 

y; / ^(Pi)d W (4) 

among all m-tuples (Vi, V2, V m ) of functions £ ^(fa) for which 

Ei=i Vi(Pi) < H(pi, for all ( Pl , ... jPm ) G R m . (5) 

Furthermore, the maximum value in ([4]) coincides with the minimum value in ([2j) and 

m—l 

Viipd = H(pu—,Pm) for all (pi, ...,p m ) € support^). (6) 

In a certain sense, the dual problem Q provides a decoupling of the original Monge-Kantorovich problem. In 
this note, we shall show that the above scheme also provides a decoupling of the system ([lj, at least for certain 
type of solutions. Roughly speaking, under certain monotonicity conditions on a solution u — (ux, u 2 , u m ) 
of ([1} , the duality in the Monge-Kantorovich problem applied to a suitable set of marginals (/x Ul , /i„ 2 , . . . , fj, U2 ) 
associated to u, leads to a decoupled system 



A, 



Ui = -^(ui(x)) for i = l,...,m. (7) 

having u — (u-i, U2, u m ) as a solution, where Vi,...,V m are the corresponding Kantorovich potentials. 
Decoupled systems are much simpler than coupled systems and have a number of advantages. For example, 
whenever the decoupling above is possible and the potentials Vi are positive, one can deduce the following 
Modica inequality [10] for systems 

E™ilV^(^)| 2 <2H( Ul (x),...,u m (x)) for all x e R N . (8) 

Whether such an inequality holds true for a general gradient system remains an open problem. See Caffarelli- 
Lin [4] and Alikakos [I]. 

The compatibility condition is also equivalent -up to a change of variables- to the classical notion of 
submodularity (also known as 2-monotonicity) of H (see Definition ^. 3p . A result of Carlier on the structure 
of optimizers in ([2]) for a submodular cost function H is essentially equivalent to a rearrangement inequality 
of Hardy-Littlewood type [S] [5] [3] • We shall see that the covariant analogue of this result noted in [TT] implies 
a rearrangement inequality for compatible H, but where decreasing rearrangements should be replaced by 
H -monotone rearrangements. Note that rearrangement inequalities have found applications in nonlinear 
optics 8, , where one looks for ground states of energy functionals of the form 

r i m 

E(u)= / -V IVu^z)! 2 +H(u)dx. (9) 
it" 2 i=1 

Note that ([T]) is nothing but the Euler-Lagrange system corresponding to this energy functional. 
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For a submodular H, any nonnegative function u may be rearranged in a symmetric and decreasing way 
(that is, the decreasing rearrangement is with respect to the variable |x|) to decrease the total energy. In our 
case, where H is orientable, any u can be rearranged in an H- monotone way with respect to the variable xjv 
to decrease the second component of the energy. In N — 1 dimensions, this rearrangement will also decrease 
the gradient term, so that the total energy also decreases under iJ-monotone rearrangement. 

Finally, we note that, while these types of decouplings can only be done for gradient systems, it is not 
essential to have the Laplacian on the left hand side; similar decouplings are possible when the left hand 
side is replaced by any decoupled differential operator DiUi. 

2 Orientable systems, compatible cost and sub-modular functions 

We begin by recalling the definitions of an orientable system and of a compatible cost. Our definitions here 
are actually strict versions of the original definitions in [BJ, where non-strict inequalities are used. 

Definition 2.1. The system (or the non-linearity H) is called orientable if there exist constants B{ such 
that, for all i ^ j, 

„ „ d 2 H 

<0. 



dpidpj 

Next, we recall the definition of compatibility, discussed in |11) . 
Definition 2.2. We say H is compatible if for all distinct k, we have 

& < °f° r ««P= (Pi.Pa.--.Pm) • 
Next, recall the following classical definition. 
Definition 2.3. The function H is submodular (or 2-increasing in Economics) if 

H(x+ he l + kej) + H(x) < H(x + he l ) + H(x + kej) (i ^ j, h,k>0), 

where x = (x±, ...,x m ) and denotes the i-th standard basis vector in W n . 
Lemma 2.4. The following are equivalent for a function H £ C 2 (]R m ,R). 

1. H is orientable. 

2. H is compatible. 

3. After a change of variables, H is submodular. 

Proof. We will first show that 1) and 2) are equivalent. First suppose H is orientable. Then we have, for all 
distinct i,j,k, 



and therefore 



dpidpj \dpkdpj > dpkdpt 
n0 d 2 H ( d 2 H \-i d 2 H 



<0, 



/ crti x- 1 crn <Q 
dpidpj \ dpkdpj J dpkdpi 



which implies compatibility. 

On the other hand, if H is compatible, set 9i — 1 if ^ p .Q pi > and 9, ■■ — — I if d p.g pi < 0. Then, up to 
positive multiplicative constants, we have 

„ „ d 2 H d 2 H d 2 H d 2 H 

< 0, 



3 dpidpj dpidpi dpidp dp l dp i 
establishing i7-orientability. 
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Now we will prove the equivalence of 2) and 3). First, we recall that a smooth function H is submodular 

if 

< for all i j and p = (pi,p 2 , •■•,Pm)- (10) 

As was noted in pj], compatibility is equivalent to the existence of changes of variables pi n- qi such that 
H(qi, q2, q m ) satisfies (|10| . It is therefore invariant under this sort of transformation. Assuming now that 
H is compatible, define a change of coordinates as follows: set q\ = p\ and, for i > 2, set qi = pi if q Vi q v . < 

and <?i = — pi if ^ x Q p . > 0. It is then clear that q q < for all i. For any distinct i,j ^ 1, we then have 
by the compatibility condition, 

d 2 H ( d 2 H \-i <9 2 # 

<0, 



\dq\dqj J 



dqidq l 



which easily yields Q q .g q . < 0. 

On the other hand, to see that 3) implies 2), it is sufficient to note that submodularity implies compati- 
bility, and that compatibility is invariant under changes of coordinates of the form pi t— > qi(pi). □ 

We note that the condition of submodularity on H is essentially equivalent to the following extended 
Hardy-Littlewood inequality: for all choices of real- valued non- negative measurable functions (u\, ...,u m ) 
that vanish at infinity, we have 

H(u\(x), ...,u* m (x))dx < I H(ui(x),...,u m (x)dx, (11) 

where u* is the symmetric decreasing rearrangement of Ui for i = 1, m. See for example Burchard-Hajaiej 
[3] and the references therein. 

We now describe the properties of multi-marginal mass transport under a compatible function cost H. 
Given probability measures /ii,/i2, ...,/im on R, and if g d q .g q . < (i.e., H is submodular relative to the q 
variables), then as shown by Carlier [5], there is a unique solution to the optimal transportation problem ([2]). 
given by 7 = (/, fo, fs, •••/m)#A*i ) where fa : R — > K is the unique increasing map pushing forward fi\ to /i^. If 
now H is compatible in the original p coordinates, then the unique solution to (J3]) is 7 = (Id, 51, 32, Sm)#/^i, 
where q, ; : R — > M is increasing if „ a £ < and decreasinq if a a > 0. It is also the unique such map 
pushing forward /ii to (fxi,fj,2, Mm)- For proofs and more discussion of the compatibility condition, see 

HQ 

We now discuss the notion of 7i-monotonicity introduced by Fazly-Ghoussoub in [BJ. 

Definition 2.5. 1. A function u — (u%, u m ) G C^R^R" 1 ) is said to be monotone if each Ui is 

strictly monotone with respect to xn ; that is, if J-^- ^ 0. 

2. u is said to be H -monotone if it is monotone and if for all i j , 

d 2 H dui duj 



dpidpj dxM 8xm 



< 0. 



It is easy to see that the existence of an 7i-monotone function implies that H is necessarily orientable. 
Conversely, we shall see below that if H is orientable, then one can associate to any function u — (u\, u m ), 
its i7-monotone rearrangement. But first, let's note the following easy application of the above lemma coupled 
with Carlier's result. For our purposes it will often be useful to decompose x G R N into (2/, xn) G R^ -1 x R. 

Proposition 2.6. Let u = (u\, U2, u m ) be a bounded monotone function in C 1 (R Ar ;R m ). Let fi be a 
probability measure on R that is absolutely continuous with respect to Lebesgue measure. For each x' € R^^ 1 , 
let fj,f be the pushforward of /1 by the map uf : xn i— > Ui(x',X]y) and set 7^ := (uf ,u% , it^J^/x. Then 
the following are equivalent: 

1. u is H -monotone. 

2. For each x' G R^ 1 , the measure j x is optimal for the Monge-Kantorovich problem (|2|). when the 
marginals are given by /m? . 



4 



Proof. 1) — > 2) is obvious. Assuming now 2), we note that for each x' <E i^ -1 , the measures fif are 
absolutely continuous with respect to Lebesgue measure by the (strict) monotonicity of the tij. The support 
of the optimizer j x must be _ff-monotone, by the equivalence of H orientability and H compatibility, 
combined with the result in [llj . As this support is exactly the image of (uf , uf , ...,it^J, this implies 



d 2 H dm du 



< 0. 



dpidpj dxM 8xm 

As this holds for any x' , this completes the proof. □ 

The above proposition leads to the following notion: 

Definition 2.7. For a function u — {u\,U2, ...,u m ) in C(WL N ;M. m ), we define its H -monotone rearrangement 
u*(x' 7 xn) = (u*, u%, u* n ) with respect to x' as follows: For each fixed x' £ R* -1 , denote by /if the 
push-forward of Lebesgue measure /i on M by xn <— > Ui(x' 'xn), and define 

• For i = 1, xjv | — ^ u[(x' , x n) to be the unique monotone decreasing map which maps [i to \i\ . 

• For i 1, then 

— If Qp.g pi < 0, let x n i— > u*(x',xn) be the unique monotone increasing map pushing /j, to fj,f . 

— If Qp.Q pi > 0, let xn u*(x' , xpf) be the unique monotone decreasing map pushing /i to /if . 

By combining the above Lemma with a change of variables and the rearrangement inequality for 
submodular functions, we obtain that if H is orientable, then 



H (u*(x', xjsj))dxM < j H(u(x', xn))oIxn ■ 

Jk 

Integrating over x 1 then yields 

H{u* {x' ,xjsi))dx < I H (u(x' , xw))dx. 



That is, the _ff-monotone rearrangement lowers the second term in the energy functional E(u) in equation ([9]). 
As to the gradient term in E(u), it is known that it decreases under monotone symmetric rearrangements. 
On the other hand, it does not necessarily under H- monotone rearrangements, unless N = 1. We therefore 
have the following. 

Proposition 2.8. If H is orientable, then for any u = (ui, tt2, u m ), we have 

1. J RN H(u*(x))dx < J RN H(u(x)dx, where u* is the H -monotone rearrangement of u. 

2. If u is one- dimensional, i.e. u £ C 1 (R;R m ), then E(u*) < E(u), where 



E <»> = / R (5 g i^l" 



and consequently, the system of ODEs, 

u xx = VH(u) on R, (12) 
has an H-monotone solution, whenever E attains its infimum on _ff 1 (R;R m ). 
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3 Decoupling systems in the presence of i7-monotone solutions 



Fazly and Ghoussoub [6] conjectured, and proved in dimensions N < 3, that if-monotone solutions of the 
system of elliptic PDEs on R N , 

Au = VH(u), (13) 

where u — (ui, u 2 , u m ) represents an m-tuple of functions, must be essentially 1-dimensional; that is, each 
Ui takes the form Ui(x) = gi(cti ■ x' — x^) for some € IR^ -1 . This is a systems analogue of DeGiorgi's 
famous conjecture for monotone solutions of the Allen-Cahn equation [7] [2]. 

They also showed that, for orientable systems in dimension N = 2, all the components of a stable solution 
u = (u\, U2, u m ) have common level sets, which also happened to be hyperplanes. We shall therefore say 
that the components («i)i=i have common level sets if for any i ^ j and any A g R, there exists some A such 
that 

{x : Uiix) = A} = {x : Uj(x) = A}. 

One might therefore conjecture that stable solutions of orientable systems may always have common level 
sets in higher dimension N, if not necessarily hyperplanes. This can be seen as complementary to the 
conjecture of Fazly and Ghoussoub, which asserts that -at least in low dimensions (say N < 8) - the level 
sets of iJ-monotone solutions are hyperplanes (but possibly different hyperplanes for each component Ui of 
the solution). 

We now establish a result on the decoupling of the system (fT3|) . 



Theorem 3.1. Let u = (ui, ...,u m ) be a bounded solution to system (|13p . 

1. If u = (iti, ...,u m ) is monotone, then there exist functions Vi(pi,x') such that u solves the system of 
decoupled equations: 

Au i = f^Mz),0 fori=l,...,m. (14) 
Furthermore, along the solution, we have 

YaLi Vi(ui(x),x') = H(ui(x),u 2 (x), ...,u m (x)) for x e R N . (15) 

2. If u ~ (ui, u m ) is H -monotone, then for all p = (pi,p2, ...,p m ) € R m , 

m 

,P2,-,Pm)- (16) 

3. If the Ui have common level sets, then the Vi can be chosen to be independent of x' , that is Vi(pi, x') = 
Vi{ Pi ). 

Proof. Fix x' 6 R^ -1 , and define Vi(-,x') on the range of xn > Ui(x' ,xn) as follows. Choose the unique 
xjv such that pi = Ui(x' ,xn) and set 

dV l , dH , 

OPi OPi 

It then follows by construction that 

. dH . dVi . 

Am = ——(u[x ,x N )) = ——(Ui,x ). 
opt dpi 

Note that, as each Vi(x' , ■) is defined only up to an arbitrary constant, we can choose the constants so that 
for xn = 0, 

m 

2J Vi(ui(x', 0),x') - H(ui{x' , 0), u 2 (x' , 0), u m (x', 0)) = 0. 
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Now note that by the definition of Vi , we have 
d 



} J Vi(Ui(x', x N ), x') - H{u\{x , x N ), u 2 {x' , xn), —, u m {x , x N )) 



8xn 

EdVi , , dui ^ dH , , , dui 
— {Ui{x ,X N ),X )- } ——(Ui(X ,XN),U2(X , X N ), ... , U m [X ,X N ))- 
, =1 dpi dx N ^ dpi dx N 

= 0. 

Therefore, we have 

m 

^ Vi(ui(x' ,x N ),x') - H(ui(x', xn), u 2 (x',xn), ■ u m (x' ,xn)) = 

everywhere on M. N . 

2) Now suppose that u is ii-monotone. For fixed x' S R^ -1 , let 7 X and {/if ,i = l,...,m} be defined 
as in Proposition 12.61 The measure "f x is then optimal for the Monge-Kantorovich problem ([2]) with given 
marginals /if , i = 1, m. In this case, the Vi defined above play the role of Kantorovich potentials, and so 
we have 

m 

22Vi{pi,x') < H(p 1 ,p 2 ,...,p m ). 

i=l 

everywhere. 

3) In this case, the image of the map (uf , uf , u^) is independent of x' . The measures j x are then all 
supported on the same set, and so we can choose the Vi(pi,x') = Vi(pi) to be independent of x' . The image 
of (ui, U2, u m ) is equal to the image of (uf ,u x , u^) for any fixed x' , and the results in jTT] imply that 
any measure supported on this set is optimal for its marginals. □ 

Remark 3.2. If u is i?-monotone such that all its components (ui)i have common level sets, then the 
image of (ui,U2, ...,u m ) is a 1-dimensional set, and, by iJ-monotonicity, the results in [11] imply that any 
measure supported on this set is optimal for its marginals. In particular, for any measure n on K , absolutely 
continuous with respect to Lebesgue measure, the measure 7 = (ui, u 2 , w m )#M is optimal for its marginals 
l-h = Ui#(i. 

As an immediate application of Theorem 13. 11 we obtain the following Modica type estimate. 

Corollary 3.3. Suppose that u is a bounded monotone solution of (j 1 3[) such that all its components (ui)i 
have common level sets. Then, letting Vi = min p . gflarase ( u .) Vi{pi), we have 

m rn 

-£|v Ui (*oi 2 <ffM-£^ (17) 
i=i 1=1 

Proof. By Theorem 13. 11 u solves the decoupled system 

A«i = ——(Ui). 
dpi 

Then, by Modica's well known inequality, we have ^|Vui(a:i)| 2 < l^(uj) — V,. Summing over i and using 
Theorem 13. 11 part 1), completes the proof. □ 

Remark 3.4. This corollary applies to, for example, bounded, monotone solutions in IR 2 with an orientable 
H. Fazly and Ghoussoub [5] showed that in this case bounded stable solutions satisfy 

Vitj = QjVui (18) 



See [6] for the definition of stability. 
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for constants Cij with signs opposite the signs of d p.Q p . ■ When the level sets of Ui (x', x n) — Qiifli'd —%n) 
are hyperplanes, this implies that dj = a,j := a for all Therefore, the Mi's share the same level sets and 
the corollary applies. 

In fact, in this case we can prove even more; we can remove the V, terms from the inequality. A simple 
computation yields 

(|a| 2 + l)g'{(a ■ x' - x N ) = Au,(x) = V t '(gi(a ■ x' - x N )). 
Denoting the argument a ■ x' — xn by y, multiplying the equation by g'(y) and integrating, we obtain 

\{\a\ i + l)g' i {yf = V i {g i {y))-C i , (19) 

for some constants Cj. As the solution gi(a ■ x' — xn) is bounded and monotone, we can take the limit as 
y — > oo; as we must have lim^oo g'(y) = 0, we obtain Cj = V^lim^oo gi{y)). Now, summing over i in (|19|) 
yields 

1 m m m m 

-(|a| 2 + i)^^(y) 2 = E^^(y))-E^ = i/ (-9i(y)'52(y),...,5 m (y))-E c '- ( 20 ) 

i—l i—1 i—l i—1 

Taking the limit as y oo of the equality 

rn 

H(gx(y),g 2 (y), -,g m (y)) = ^2v{gi(y)) 

i=l 

yields 

m 

H( lim lim g 2 (y),-, hm 5m(y)) = 

y—>oo y-^oo y— >oo * — » 

If is nonnegative, (|20)) yields 



i=l i=l 

< H(gx(y),g 2 (y),...,g m (y)) 
= H(ui(x),U2(x),...,u m (x)). 

Remark 3.5. It is not clear to us whether the approach above provides any more information about 
.ff-monotone solutions when the solutions are not 1-dimensional, but surprising it does imply more about 
solutions with the opposite geometry. Suppose that — H is orientable, and that u is a bounded — H monotone 
solution to (|T3"|) ; in other words, for all i ^ j 

d 2 H dui duj 

> 0. 



dpidpj 8xn dxN 



Then for each x' , the measure ^ x maximizes (rather than minimizes) the multi- marginal Kantorovich func- 
tional ([2]) and the functions Vi from Theorem 13.11 satisfy the constraint 



Vi(x', P i) > H(pi,p 2 ,...,p m ). (21) 



If the Ui also have common level sets, then the Vi are independent of x\ by Theorem 13. II If H is everywhere 
nonnegative, taking infimums in (|2ip yields y'., Vi > 0; combining this with the previous corollary implies 
the following Modica estimate for these types of solutions: 

1 - 

-Eiv^)i 2 <#w- 

i=l 
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